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Abstract 

We have studied dynamical responses and the Jarzynski equaUty (JE) of classical open systems 
described by the generalized Caldeira-Leggett model with the nonlocal system-bath coupling. In 
the derived non-Markovian Langevin equation, the nonlinear nonlocal dissipative term and state- 
dependent diffusion term yielding multiplicative colored noise satisfy the fluctuation-dissipation 
relation. Simulation results for harmonic oscillator systems have shown the following: (a) averaged 
responses of the system (x) to applied sinusoidal and step forces significantly depend on model 
parameters of magnitudes of additive and multiplicative noises and the relaxation time of colored 
noise, although stationary marginal probability distribution functions are independent of them, 
(b) a combined effect of nonlinear dissipation and multiplicative colored noise induces enhanced 
fluctuations {[x — (x)]^) for an applied sinusoidal force, and (c) the JE holds for an applied ramp 
force independently of the model parameters with a work distribution function which is (symmetric) 
Gaussian and asymmetric non-Gaussian for additive and multiplicative noises, respectively. It has 
been shown that the non-Markovian Langevin equation in the local and over-damped limits is 
quite different from the widely adopted phenomenological Markovian Langevin equation subjected 
to multiplicative noise. 
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I. INTRODUCTION 



In the last almost half a century, many studies have been made on the Langevin model 
which is widely employed as a useful model for a wide range of stochastic phenomena (for 
a recent review, see Ref. [l]). Dynamics of a Brownian particle subjected to potential V{x) 
is modeled by the Langevin equation given by 

X = -V'{x)-^oX + V2Dat), (1) 

where dot and prime stand for derivatives with respect to time and argument, respectively, 
7o denotes dissipation, C,(t) is zero-mean Gaussian white noise with (^(t)^(t')) = — and 
D expresses the strength of noise. Dissipation and diffusion terms satisfy the fluctuation- 
dissipation relation (FDR), 

D = ksT-fo, (2) 

where ks is the Boltzmann constant and T the temperature. The FDR implies that dissipa- 
tion and diffusion processes originate from the same event. The simple Langevin equation 
given by Eq. ([T]) is based on the two assumptions: (a) a dissipation is local in time and (b) 
a diffusion depends on velocity x(t) but is independent of state x(t). State-independent and 
state-dependent diffusions are commonly referred to as additive and multiplicative noises, 
respectively. Multiplicative noise can be phenomenologically described in a number of ways: 
for example, a diffusion term in Eq. ([1]) may be generalized as (for a review of study on 
multiplicative noise, see Ref.[2|) 

V2Dm ^ V2DG{x)m, (3) 
^ V2A^^it) + V2MGix)Ut), (4) 

where G{x) is a function of x, A and M denote strengths of additive and multiplicative noises, 
respectively, and ^i(t) and ^2(^) are white noises. However, phenomenological diffusion terms 
given by Eqs. ([3]) and (jlj) have no microscopic bases and the FDR for such diffusions is not 
definite |3i]. The stationary probability distribution function (PDF), which is obtained from 
the Fokker-Planck equation (FPE) corresponding to the Langevin model including these 
diffusion terms, is generally different from the Boltzmann factor, Q-^m/i^BT 

The importance of going beyond the assumptions (a) and (b) has been recognized over 
many decades. The microscopic origin of additive noise has been proposed within the frame- 



work of system-bath Hamiltonians [7H9|, which are known as Caldeira-Leggett (CL) type 
models. By using the generalized CL model including a nonlinear system-bath coupling, we 
may obtain the non-Markovian Langevin equation with nonlocal dissipation and multiplica- 



tive noise which preserves the FDR 
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The nature of nonlinear dissipation and multi- 
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15| . Ref. [l2| discusses 



plicative noise has been recently explored with renewed interest 
a possibility of observing a quantum current in a system with quantum state-dependent dif- 
fusion and multiplicative noise. Dynamics in a metastable state which is nonlinearly coupled 
to bath driven by external noise has been studied Ref. [l^ investigates a temporal 

development in an average velocity of noninteracting Brownian particles in a finite system 
with nonlinear dissipative force. Quite recently, a detailed comparison is made between non- 
Markovian and Markovian Langevin equations including additive and multiplicative noises 



15| . It has been shown that in many cases, the Markovian (loc al) approximation is not a 



reliable description of the non-Markovian (non-local) dynamics [15^]. Nonlinear dissipation 
and multiplicative noise have been recognized as important ingred ients in several fields such 



as mesoscopic scale systems 16|, ll7|] and ratchet problems 181-121 1. 



In the last decade, we have significant progress in experimental and theoretical under- 
standing of nonequilibrium statistics (for reviews, see Refs. 22|-|24|). At the moment we 



have three kinds of fluctuation theorems: the Jarzynski equality 



transient-fluctuation theorem 



2J28 



E) |25j . the steady- and 



I , and the Crooks theorem [27|, |28| . These theorems are 
applicable to nonequilibrium systems driven arbitrarily far from the equilibrium state. In 
this paper, we pay our attention to the JE which was originally proposed for a classical 
isolated system and open system weakly coupled to baths [25|, |29[. Subsequently Jarzynski 
proved that the JE is valid for strongly coupled classical open systems [SOJ. A validity of 
the JE has been confirmed by various experiments for systems which may be described by 



damped harmonic oscillator models |31l-l36|. Stimulated by these experiments, many the- 
oretical analyses have been made for harmonic oscillators with the use of the Markovian 
Langevin model 33l-l36[|. the non-Markovian Langevin model 37H40l|. Fokker-Planck equa- 



tion |41|, and Hamiltonian model |42l-l47l|. The validity of the JE has been^xamined for 



anharmonic oscillators 
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48| and for van der Pol and Rayleigh oscillators [49]. We should 



note that these studies have been made for the non-Markovian and/or Markovian Langevin 
models with additive noise. Recently the JE in the Markovian Langevin model with mul- 
tiplicative white noise for Brownian particles has been discussed in Ref. js^. However, a 
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study on the JE for the non-Markovian Langevin model with multiphcative colored noise is 



scanty at the moment 



5l|. 



The purpose of the present paper is twofold: (1) to make a detailed study of the non- 
Markovian Langevin model derived from the generalized CL model for classical open systems 
with nonlinear nonlocal dissipation and state-dependent diffusion and (2) to calculate re- 
sponses to applied forces and examine a validity of the JE in the system. In the following 
Sec. II, we derive the non-Markovian Langevin equation, adopting the generalized CL model 
including nonlinear system-bath coupling IJQ (Sec. IIA). The Ornstein-Uhlenbeck (OU) 
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52| 



m 



process of colored noise is taken into account. By using the two methods 
which new variables are introduced into the non-Markovian Langevin equation, we obtain 
a set of four first-order differential equations and the relevant multi-variate FPE. The local 
limit of the non-Markovian Langevin equation is examined (Sec. IIB). In Sec. Ill, we study 
harmonic oscillator systems, applying simulation method to the four differential equations 
mentioned above. We calculate the stationary marginal PDF of the system (Sec. IIIA) 
and its responses to applied sinusoidal and step forces (Sec. IIIB). In particular, frequency- 
dependent responses of the mean position of {x(t)) to sinusoidal force have been made in 
detail. We obtain enhanced fluctuations of {[x{t) — (x(t))]^) induced by a combined effect 
of nonlinear dissipation and multiplicative noise. Applying a ramp force to the system, the 
validity of the JE has been examined (Sec. IIIC). In Sec. IV we discuss the over-damped 
limit of the Markovian Langevin model subjected to multiplicative noise. Sec. V is devoted 
to our conclusion. 



II. THE SYSTEM-BATH MODEL 



A. Non-Markovian Langevin equation 



We consider a system of a classical oscillator coupled to a bath consisting of A^-body 
uncoupled oscillators described by the CL model jo], Q , 

H = Hs + Hb + Hi, (5) 
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with 



Hs 

Hb + Hj 



p 

N 

E 

n=l 



+ V{x)-xf{t), 



2m„ 2 



(In 



(6) 
(7) 



Here Hs-, Hb and if/ express Hamiltonians of the system, bath and interaction, respectively: 
X, p and V{x) denote position, momentum and potential, respectively, of the system: 
m„ and ujn stand for position, momentum, mass and frequency, respectively, of bath: the 
system couples to the bath nonhnearly through a function 0(x): f{t) expresses an applied 
external force. The original CL model adopts a linear system-bath coupling with (j){x) = x 
in Eq. ([7]) which yields additive noise IqI. By using the standard procedure, we obtain the 



generalized Langevin equation given by [9 

m = 



V'{x{t)) - <p'{x{t)) I 7(t - t') <P'{x{t')) x{t') dt' + ct>'{x{t)) at) + fit), (8) 





with 



lit) 

at) 



N 

E 

n=l 

N 

E 

n=l 



COS LJnt^ 



g„(0) - 0(a;(O)) 



(9) 



cos u„t 



CnPniO) 



rrinUJr, 



sinUnt } , (10) 



where 7(t — t') denotes the non-local kernel and (it) stands for noise. Dissipation and 
diffusion terms given by Eqs. ([9]) and ( fTOj) . respectively, satisfy the second- kind FDR, 



(CWC(O)o = kBT^it-t'), 



(11) 



where the bracket {■)q stands for the average over initial states of q'„(0) and Pn(0) 10l-ll2|. 
We have adopted the OU process for the kernel 7(t — t') given by 



Jit - t') = e-(*-*')/^ 



(12) 



where r and 70 stand for the relaxation time and strength, respectively, of colored noise. 
The OU colored noise may be generated by the differential equation, 

Cit) , V2A;bT7o 



at) 



at), 



(13) 
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where ^{t) expresses white noise with 



m) = 0, mm) = 5it-t'). 

Equations (fT3|) and (fT4|) lead to the PDF and correlation of colored noise given by 

P(C) oc e-(^-/2^°)^', 



(14) 



{cmt')) 



e-^'-^'y^ = ksT ^{t - t') 



(15) 
(16) 



where /3 = 1//cbT. 



^. The method A 



The two methods have been proposed to transform the non-Markovian Langevin equation 
given by Eq. (|8]) into multiple differential equations [ll|, [l5, [52|. In the method A, we 
introduce a new variable u{t) [l5, 52|, 

u{t) = - [ -fit- t')(j)'{x{t'))x{t') dt', (17) 



- / -fit- t')(t)'{x{t'))x{t') dt', 
Jo 

to obtain four first-order differential equations for x{t), p{t), u{t) and ({t) given by 
x{t) = p{t), 

Pit) = -V'(a;)+0'(x(t))w(t) + /(t) + 0'(x(t))C(t), 
uit) = 



-!^-(^)0'(.(t))p(t), 

gt) ^ yw^o ^^^^^^ 



(18) 
(19) 
(20) 

(21) 



T T 

From Eqs. (|T8ll -( |2T]) . we obtain the multi-variate FPE for distribution of P{x,p, u, (, t), 
dt = -^pP(x,P,«,C,t) 



d 



-V\x) + f(t) + <p\x)u + <P\x)Q P(x, p, u, C, t) 



+ 



dp 
d_ 

du 
kBT-fo\ 



^ + (7) Pi^^P^ ^' C, t) + ^ ) P{x,p, u, C, t) 



■P{x,p,u,C,t)- 



(22) 
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2. The method B 



In the method B, we introduce a new variable z{t) 11 1 



z{t) = - (^-j - t')ct>'{x{t')) dt' + [- ] at) + m, 

to obtain four first-order differential equations, 



x{t) 
Pit) 

m 
m 



Pit), 

dU{x, z) 



dx 
1 dU (x, z 



+ fit), 



7o 



dz 



+ 



at) , V2A:bT7o 



r 



7o 



with 



U{x,z) = y(a;)+(^) [z-<P{x)] 



(23) 

(24) 
(25) 

(26) 
(27) 

(28) 



It is noted that white noises in Eqs. (126|) and (l27j) come from the same origin. 

A variable ({t) in Eq. ( l27l) is isolated from the rest of variables in the four differential 
equations. From Eq. fl27p we may obtain its stationary PDF, -P(C), given by Eq. (IT^ . The 
FPE relevant to Eqs. ( l2^ -(l26l) for P{x,p, z,t) is given by 



dP{x,p, z) 
di 



•9 . d 

-^pP{x,p,z) + — 
ox op 



dU{x, z) 
dx 



-fit) 



n/ X 1 "9 dU(x, z) , 
F{x, p, z) H F{x, p, z) 



7o dz dz 



+ 



■P{x,p,z) 



7o dz"^ 

The stationary PDF of Eq. with fit) =0 is given by |ll| 

Pix,p,z) oc e-'5[^''/2+^(-.-)], 

which leads to stationary marginal PDFs, 



(29) 



(30) 



P{x) = J P{x,p,z) dpdz oce-'^^^^'l 
P[p) = [ P{x,p, z) dx dz oc e-f^p'/\ 



(31) 
(32) 

Equations ( !T8|) - (pT|) in the method A are equivalent to Eqs. (12^ - (127|) in the method B. 
An advantage of the method A is that the local limit of r — is easily obtainable in Eqs. 



(fT8|) -( l2T]) . while in the method B the analytical expression for the stationary PDF given by 
Eq. fl30|) may be derived. In our simulations to be reported in the following section, we have 
mainly employed the method A, whose results are partly checked by separate simulations 
using the method B. 

B. Markovian Langevin equation 

It is worthwhile to examine the local limit of Eq. with a kernel 7(t — t') given by 

7(t-t') = 27o5(t-t'), (33) 

which leads to the Markovian Langevin equation, 

x{t) = -V'{x{t)) - 7o0'(a:(t))' x{t) + y^2k^o <P'{x{t)) ^(t) + f{t). (34) 

It is evident that the Markovian Langevin equation becomes a good approximation of the 
non-Markovian one in the limit of r — )■ 0. 

From Eq. flMl) . we obtain three differential equations, 

x{t) = p{t), (35) 

p{t) = -V'{x) - ^o<j>'{x{t) )Mt) + fit) + <P'{x{t)) at), (36) 

C(t) = -^ + ^^?^e(t). (37) 
r T 

The PDF for ( is given by Eq. f|T5|) . The relevant FPE for the PDF of P{x,p, t) is expressed 

by 

^E^^M. = -^pP^x^p^t) + -^[V'{x)-f{t)+jo^\x)'p]P{x,p,t) 

+ fcBT7o0'(x)2^0'(x)^0'(x)P(x,p,t). (38) 
The stationary distribution of Eq. (|38l) with f{t) = is given by 

Pix,p) oc e-^[p'/2+y(x)]_ (39) 

This is consistent with the result of the non-Markovian Langevin equation given by Eq. 
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III. HARMONIC OSCILLATOR SYSTEMS 



A. Stationary marginal PDF 

Simulations have been performed for harmonic oscillator systems where V{x) and 4>{x) 
in Eqs. ©-(IT]) are given by 

Vix) = (40) 

2 

0(a;) = ^- + bx, (41) 

where Ug stands for oscillator frequency of the system, and a and b denote magnitudes of 
multiplicative and additive noises, respectively. We have solved Eqs. f fT8|) - fl2T]) by using the 



Heun method [53|] with a time step of 0.001 for parameters of Us = 1.0, 70 = 1.0 and 
1.0 otherwise noticed. Simulations have been made for < t < 1000.0, which are averaged 
over 10"^ sets of initial states of Gaussian-distributed a;(0) and p(0) with (a;(0)) = (p(0)) = 
and (p(0)^) = (u;^x(0)^) = ksT. In all simulations, we have used the initial conditions of 
u{0) = and C(0) = 0. 

First we show marginal PDFs of P{x), P{p), P{u) and -P(C) for f{i) = 0, which are 
evaluated by simulations of Eqs. f|T8|) - fl2T|) with discarding initial results at t < 200.0. 
Figures [l](a), (b), (c) and (d) show P{x), P{p), P{u) and P(C), respectively, obtained for 
(a, 6, r) = (1, 0, 10.0) (sohd curves), (1, 0, 0.1) (dashed curves), (0, 1, 10.0) (chain curves) and 
(0, 1,0.1) (dotted curves), where r stands for the relaxation time of OU colored noise. We 
note that P{x) and P{p) in Figs. [T](a) and (b) are independent of the parameters of (a, 6, r) 
although P{u) and P(C) in Figs. [I](c) and (d) depend on them. Calculated P{x) and P{p) 
are in good agreement with Gaussian PDFs given by Eqs. (I3T1) and (!32|) . Equation (fT5|) 
shows that -P(C) is the Gaussian PDF whose variance depends on r for fixed 70 and T. In 
contrast. Fig. [T](c) shows that P{u) is Gaussian PDF for additive noise but non-Gaussian 
PDF for multiplicative noise: the kurtosis of P{u) defined by 

- (42) 

is r]^ = 3.0, 3.0, 7.8 and 10.3 for (a, b, r) = (0, 1, 0.1), (0, 1, 10.0), (1, 0, 0.1) and (1, 0, 10.0), 
respectively, where the bracket (■)„ denotes an average over P{u). Note that rju = 3.0 for 
the Gaussian distribution. 
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B. Responses to applied forces 



1. Sinusoidal forces 

Dynamical responses of harmonic oscillator systems to applied sinusoidal and step forces 
are studied. We first apply a sinusoidal force given by 

f{t) = g sin (^^^ = g sin Wot, (43) 

where g, Tq and Uq (= 27r/To) denote its magnitude, period and frequency, respectively. A 
sinusoidal force given by Eq. f H5]) yields averaged outputs given by 

= (44) 

fxpit) = {pit)) = fi^it), (45) 

Fourier-transformed outputs given by 

POO 

fi4u] = / e'^'fi,it)dt, (46) 



'OO 

oo 

i uit , 



fiploo] = / e'^'fi^,{t)dt, (47) 

J — oo 

have peaks at u = Uq. Output magnitudes defined by 

4(wo) = lyUxMP, (48) 
Ipiujo) = IfipluoW^ = ujq Ixi(^o), (49) 

express frequency- dependent responses of the system. 

In the case of additive noise only [(j){x) = x], we may obtain the analytical result of Ixi^oo) 
with the susceptibility given by 

RxM = = — ^ : r^, (50) 

where 'y[u] is the Fourier transform of the kernel, 7(t), 

7M = r^^- (51) 

1 — lUT 

From Eqs. fH51) - flFIl) . the output magnitude is given by 



nVil + coy) 



[{uji-u;^,)^{l + u;iT^)+uiji] 
7rV(l + a;oV) 
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for uo > 0, (53) 
for uq = Ug. (54) 



Equation (15^ shows that Ix{<^o) at a resonance frequency {uq = Us) is monotonously in- 
creased with increasing r. 

In the case of a general yielding multiplicative noise, however, we cannot make an 
analytical study because the Fourier or Laplace transformation cannot be employed. Then 
we have to rely on numerical simulations of Eqs. fll8p -f l2ip . Figures HJ^a) and (b) show time 
courses of fix(t) for r = 0.1 and r = 10.0, respectively, with (a, 6) = (0, 1) and (1,0) when 
a sinusoidal force given by Eq. fj43|) with g = 1.0 and Tq = 10.0 is applied. The results 
of {a,b) = (1,0) and (0,1) are almost the same for r = 0.1 in Fig. 2(a). In contrast, for 
T = 10.0, /i^.(t) of (a, 6) = (1,0) is a little different from that of (a, 6) = (0, 1) in Fig. 2(b): 
an irregularity in the former is larger than that of the latter. The irregularity in fix{t) for 
(a, b) = (0, 1) is gradually reduced at a larger t (relevant result not shown). 

Although the difference between fix{t) for different r values with additive and multiplica- 
tive noises is not so clear in Fig. 2, it becomes evident in the Fourier-transformed quantity 
of or Ix{uo). Figure |3]^a) shows the coq dependence of output magnitudes of Ix{uJo) 

with (a, 6) = (0,1) (chain curve) and (1,0) (solid curve) for r = 10.0. The chain curve for 
additive noise with (a, b) = (0, 1) has a resonance peak at the frequency of a bath oscillator 
{u = Us = 1.0). The dashed curve expresses a theoretical result for additive noise calculated 
by Eq. ( 153|) . which is in good agreement with a relevant result of simulations except for 
Uq ~ Us where our simulation overestimates Ix{uq). In contrast, the solid curve in Fig. 
IHl^a) for multiplicative noise has a A-type peak at uq = Ur ^ 1.4 which is different from 
Us- Furthermore, a shape of the solid curve is rather peculiar and different from that of 
the chain curve for additive noise. Figure [3]^b) shows a similar plot of Ix{u}o) for r = 0.1 
with (a, 6) = (0,1) (chain curve) and (1,0) (solid curve). A broad peak at cuq ~ 1 ^or 
multiplicative noise is smaller than that for additive noise. 

Figures lll^a)-(d) show the Lissajous plots of fixit) versus f(t) with multiplicative noise 
of (a, b, t) = (1, 0, 10.0) for Tq = 10.0, 5.0, 4.0 and 3.0 which correspond to uq = 0.63, 1.26, 
1.57 and 2.09, respectively. We note that fixit) is in phase with f{t) for Tq = 10.0 but in 
anti-phase for Tq = 3.0, and that the transition from an in-phase to an anti-phase occurs at 
Tq ~ 4.5 {uq = 1.4 = Ur). This is consistent with a peak position of /^(wo) for multiplicative 
noise shown in Fig. |3]^a). 
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We study fluctuations of x{t) in tlie system defined by 



(55) 



It is noted that in stationary state without forces, we obtain Px{t) = Pp(t) = ksT in 
harmonic oscillator systems because marginal stationary PDFs are given by P{x) oc e~^^^^^ 
and P{p) oc e~^^^/^ regardless of values of a, b and r. It is not the case when sinusoidal 
forces are applied to the system with multiplicative noise, as will be shown in the following. 
Figures a) and (b) show time courses of Px{t) in harmonic oscillator systems for r = 0.1 
and T = 10.0, respectively, for multiplicative noise of (a, 6) = (1,0) when sinusoidal forces 
with g = 1.0 and Tq = 6.0, 10.0 and 100.0 are applied. All ps{t) start from Ps(0) = 1.0 at 
t = 0.0, results for Tq = 6.0 and 10.0 being shifted by four and two, respectively, for a clarity 
of the figure. For To = 100.0, we obtain Px{t) — 1.0 for r = 0.1 and 10.0. In contrast, Px{t) 
for Tq = 10.0 and 6.0 is much increased than unity. This is more clearly seen in Figs. El^a) 
and (b), where the stationary value of pxs defined by 



is plotted as a function of Uq for additive (dashed curves) and multiplicative noise (solid 
curves). The solid curve in Fig. [6](a) expressing pxs for r = 10.0 with multiplicative noise 
has a peak at uq ~ Ur where Ur — 1.4 discussed before. The magnitude of a peak for r = 0.1 
with multiplicative noise in Fig. [6](b) is less significant than that for r = 10.0 in Fig. [6](a). 
On the other hand, dashed curves in Figs. [6](a) and[6t^b) expressing Px{t) for additive noise 
of (a, b) = (0, 1) are given by Px(t) = ksT = 1.0 independently of Tq and r. 

Solid and chain curves in Fig. [7] express the r dependence of pxs with Tq = 5.0 and 
To = 10.0, respectively, for multiplicative noise. With increasing r, pxs with Tq = 5.0 is 
much increased than that with Tq = 10.0 because the former is closer to 2'K/ujr ~ 4.19 than 
the latter. For additive noise, we obtain pxs = 1.0 (dashed curve) as mentioned above. These 
enhanced fluctuations arise from a combined effect of nonlinear dissipation and multiplicative 
colored noise. 



Pxs — 



Pxit) 



at t ~ 1000.0, 



(56) 
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2. Step forces 



Next we apply a step force given by 



for t < t 



1? 



(57) 



m 



g for t > ti, 



where ti is the starting time of an force with a magnitude of g. Figure [8]^a) shows time 
courses of fix{t) with additive noise with (a, 6) = (0, 1) for r = 0.1, 1.0 and 10.0 when step 
forces with g = 1.0 are apphed at t = ti = 100.0. For r = 10.0, an oscillation induced by 
a step force applied at t = 100.0 remains for a fairly long period after ti, whereas those for 
r = 1.0 and 0.1 are quickly damped out. Figure [8]^b) shows similar time courses of fix{t) for 
multiplicative noise with (a, 6) = (1,0) for r = 0.1, 1.0 and 10.0. Comparing Fig. [8]^a) with 
Fig. Et^b), we notice that an oscillation with multiplicative noise is more quickly disappear 
than that with additive noise. 

Figure E] shows /i^(t) for four sets of (a, 6) = (0.0,1.0), (0.2, 0.8), (0.5, 0.5) and (1.0, 
0.0) with r = 10.0. With increasing a component of the multiplicative noise, the oscillation 
induced by an applied step force is rapidly decayed. 

C. The Jarzynski equality 



where W stands for a work made in a system when its parameter is changed, the bracket 
means the average over the work distribution function (WDF), P{W), of a work performed 
by a prescribed protocol, and AF denotes the free-energy difference between the initial and 
final equilibrium states [Eqs. (jHS]) and fIBB]) ]. Equation flSHl) includes the second law of the 
thermodynamics, {W)w > AF, where the equality holds for the reversible process. The JE 
in Eq. (158!) may be rewritten as 



The JE is expressed by 





(58) 



w — 



AF. 



(59) 



When the WDF is Gaussian given by 



P{W) 



1 



(60) 



e 



W 
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we obtain 



R 



(61) 



with 



w 



((w^-M'> 



w 



(62) 
(63) 



where fiw and o"^ express mean and variance, respectively, of the WDF. Equation fl6Tl) is 
not vahd when the WDF is non-Gaussian. 
We apply a ramp force given by 



fit) 







for t < 0, 



gi^) for < t < Tf, 
g for t> Tf, 



(64) 



where r/ stands for a duration of the applied force with a magnitude g (= 2.0). 

The free energy difference AF between equilibrium states with f{t) = and f{t) = g {g, 
constant) is given by 



30|,|49| 
AF 



F{g)-F{0) = --\n(^ 



(65) 



with 



Zsig) 



Tr {e 



} 



(66) 



Tr{e-^^s} 

where Hs{g) denotes the system Hamiltonian with f{t) = g. After some manipulations 
(detail being given in the Appendix), we obtain 

^2 



AF 



9 



2oj 



2' 



(67) 



independently of a, h and r, which becomes AF = —2.0 for ujs = 1.0 and g = 2.0. 

Simulations have been performed with the same parameters as in Sees. Ill A and III B, 
but over 10^ sets of initial states. Simulation results are presented in Figs. [TUlfT^ Figure 
[TU T a). ITUT b) and ITUT c) show WDFs for ramp forces of Tf = 100.0, 10.0 and 1.0, respectively, 
with four sets of parameters of {a,b,T) = (1,0,10.0) (sohd curves), (1, 0, 0.1) (dashed 
curves), (0, 1, 10.0) (chain curves), and (0, 1, 0.1) (dotted curves). We note in Fig. [TOT a) 
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that all WDFs for r = 100.0 locate at jjw ~ —2.0 with widths of aw ~ 0.3 and that 
WDFs for multiplicative noise are asjTumetric non-Gaussian, while those for additive noise 
are (symmetric) Gaussian. Indeed, the kurtosis of the WDF given by 



- 7;;tt^ ' ^^^) 



is T]w = 3.0, 3.0, 3.9 and 4.2 for (a, b, r) = (0, 1, 0.1), (0, 1, 10.0), (1, 0, 0.1), and (1, 0, 10.0), 
respectively, with r/ = 100.0. When Tf is reduced to 10.0, behaviors of P(W) are changed. 
Figure [TUT b) shows that WDFs for r = 10.0 with different sets of (a, b) show different 
behavior but with almost the same values of fiw ~ —1.8 and aw ~ 1.0. WDFs for multi- 
plicative noise much departs from Gaussian distribution: the kurtosis of P{W) for r = 10.0 
is riw = 3.0, 3.0, 7.0 and 5.6 for (a, b, r) = (0, 1, 0.1), (0, 1, 10.0), (1, 0, 0.1), and (1, 0, 10.0), 
respectively. When Tf is further reduced to 1.0, we note in Fig. [rOT c) that all WDFs become 
almost identical Gaussian {rjw — 3.0) with fiw ~ 0.0 and aw ~ 2.0. 

These changes in fiw, <^w and r]w as a function of r/ are shown in Figs. ITTTa). fTTTb) and 
[TlT c). respectively, where marks express simulation results and curves are plotted only for a 
guide of the eye. With decreasing r from 100.0 to 0.1, /i^y changes from [iw ~ —2.0 to 0.0 
while aw increases from 0.3 to 2.0. The kurtosis riw of PiyV^ for multiplicative noise has 
a maximum around r ~ 10.0 whereas that for additive noise keeps ?7 = 3.0 independently 
of r. We should note that [iw and aw may show oscillations at r > 10 if simulations are 
performed with finer meshes (see Fig. 3(a) and (b) in Ref. 47(]). The Tj dependence of R 
calculated by Eqs. (!58|) and( l59|) is shown in Fig. [TTT d) where the equality: R = AF holds 
within conceivable numerical errors. The JE is expected to hold in our system independently 
of the parameters of a, b, r and Tf. 

Figure [12] shows the temperature-dependent P{W) for ksT = 1.0 (solid curve), 5.0 
(dashed curve) and 10.0 (chain curve) when a ramp force of Tf = 10.0 is applied to a system 
with multiplicative noise of (a, b, r) = (1, 0, 10.0). With increasing the temperature, a width 
of WDF is increased and its departure from the Gaussian distribution is more significant: 
the kurtosis of P(W) is rjw = 5.6, 15.0 and 13.0 for ksT = 1.0, 5.0 and 10.0, respectively. 
The shape of non-Gaussian WDF with multiplicative noise considerably depends on the 
temperature. 
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IV. DISCUSSION 



The over-damped limit of the Markovian Langevin equation is conventionally derived 
with setting x = in Eq. flM]) . This is, however, not the case where dissipation and 
diffusion constants are state dependent as in our case. Sancho, San Miguel and Diirr ^ 
have developed an adiabatic elimination procedure to obtain an exact Langevin and FPEs 
in such In order to adopt their method [4], we rewrite Eq. fl34l) as 



x{t) = -rix{t))-X{x{t))xit)+g{x{t))at) + f{t), 



(69) 



with 



X{x) = 7o0'(x)^ (70) 
gix) = v/2A;bT7o 0'(x). (71) 

By the adiabatic elimination in Eq. f l69|) after Ref. j4|, the FPE in the Stratonovich inter- 



pretation is given by 

dP{x,t) 



d 



V'ix) - fit) + ksT 



d_ 

dx 



P{x,t), 



(72) 



dt dx X{x) 

where we employ the relation: (7(x)^ = 2kBTX{x) derived from Eqs. ( ITOl) and (ITTjl . The 
corresponding Langevin equation is given by {4] 

[V'ix) - fit)] 1 



X 



g'ix)gix) + ^at). 



(73) 



A(x) 2A(x)2 ^ ^ ^ ' A(x) 

Note that the second term of Eq. (1731) does not appear when we obtain the over-damped 
Langevin equation by simply setting x = in Eq. fl69l) . It is easy to see that the stationary 
distribution Psix) of Eq. (1721) with fit) = is given by 



PJx) oc e-^^("). 



(74) 



In the case of <^(x) = ax'^/2 + bx and fit) = 0, Eqs. ([701), (HI]) and ([75]) lead to the 
Langevin equation with additive and multiplicative noises given by 

V'ix) keTa , / 2A;bT 



X 



+ 



7o(6 + ax)2 7o(6 + ax)3 y7o(6 + ax)2 
In the case of additive noise only (a = 0), Eq. (1751) becomes 

V'ix) , l2kBT 



(75) 



X 



7o62 



+ 



7o62 



(76) 
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In the opposite case of multiplicative nose only (6 = 0), we obtain 



X 



(77) 



Equation 
given by jsl, Is 



75 



is quite different from a widely-adopted phenomenological Langevin model 



X 



V'{x) + V2Ai{t) + V2Mx7]{t), 



(78) 



where A and M stand for magnitudes of additive and multiplicative noises, respectively, 
and ^{t) and ri{t) express zero-mean white noise with unit variance. The stationary PDF 
obtained from the PFE for Eq. (175]) in the Stratonovich sense is given by 

V'(x) 



\nP(x) 



dx 



\n{A + Mx^ 



(79) 



{A + Mx2) 

Equation (I79ll yields Gaussian or non-Gaussian PDF for V'{x) = x, depending on A and 
M |6[. The Langevin model given by Eq. (!78|) is one of origins leading to Tsallis's 



nonextensive statistics 



55|. 



V. CONCLUSION 



Dynamical responses and the JE of classical open systems have been studied with the 
use of the generalized CL model yielding the non-Markovian Langevin equation in which 
nonlinear dissipation term and state-dependent diffusion term satisfy the FDR [Eq. flTT])]. 
Simulation results for harmonic oscillator systems are summarized as follows: 

(i) marginal stationary PDFs for x and p are given by P{x) oc e~^^*-^^ and P{p) oc e"^^^/^, 
respectively, independently of a, b and r (Fig. [1]), 

(ii) Ix{ooo) with multiplicative noise for an applied sinusoidal force has a peculiar A-type 
peak at Uq = Ur (> u)s) while that with additive noise has an almost symmetric resonance 
peak at uq = cOs (Fig. |3]), 

(iii) enhanced fluctuations of Px{t) may be induced by applied sinusoidal force with mul- 
tiplicative noise (Fig. 5), and the uq dependence of its stationary p^s with multiplicative 
noise has a larger A-type peak at Uq ~ Ur for a larger r while p^s with additive noise is 
independent of uq (and r) (Figs. [6] and [7]), 
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(iv) dynamical responses to applied step forces considerably depend on a, b and r (Figs. [8] 
and [9]), and 

(v) the JE is valid independently of the parameters of a, b, t and Tf with the WDF which is 
Gaussian and asymmetric non-Gaussian for additive and multiplicative noises, respectively 
(Figs. [Hand eg). 

The items (i)-(v) imply that a nonlinear dissipative term and a state-dependent diffusion 
term have significant effects on nonequilibrium properties, although they have no effects on 



stationary marginal PDFs. The item (i) is in consistent with Refs. 25|, |30|, numerically 



confirming the analytical result of Ref. ll|. The item (i) is, however, in contrast to the 
result of a widely adopted phenomenological Langevin model subjected to additive and 
multiplicative noises [e.g., Eq. fl75]) ]. in which stationary PDFs depend on magnitudes of 
the two noises. An asymmetric non-Gaussian WDF in the item (v) is similar to that in 



the system with anharmonic potential where the JE holds [38[. Theoretical results in this 
paper are expected to bear valuable consequences on the field of open systems with non- 
local dissipation and state-dependent diffusion. It would be possible to apply the present 
approach to classical open systems such as free particles, anharmonic oscillators and bistable 
ones, and to extend it to quantum open systems. 
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Appendix: A. Free energy difference of the system 

We will calculate the free energy difference of AF, evaluating the system partition func- 
tion Zs{g) given by Eq. fl66|l . 

Zsi9) = (Al) 



where 



= Tr{e-'^^-}, (A2) 
Z{g) = Tr{e-'^[^s(»)+^^+^^]}, (A3) 
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with 



H 



B 



N 

E 

ra=l 



2m„ 2 



2 2 
Hb + Hj = Y, 



n=l 



2m„ 2 \ ninUJn 



(A4) 

(A5) 
(A6) 

(A7) 



Hs{g) denoting the system Hamihonian given by Eq. ([6]) with f{t) = g. By using the 
Gaussian integral, we obtain the classical partition function Zb in Eq. (IA2p for A^-body 
uncoupled harmonic oscillators given by 

^ POO 



N 

n 

n=l 



n=l 



n I ""in 
00 J —00 



-/3(p2/2m„+m„a;2g2/2) 



27r 

/3Un 



(A8) 
(A9) 



where h is the Plank constant [56]. Similarly Z{g) in Eq. ( 1A3I) is given by 



00 ■/ —00 

00 /"OO 



n=l 



(AlO) 



N 



dp / dy e 



-00 J —00 

-'^ f>00 /"OO 

XII/ dVr, / rfWr, e-^(P'/2m„+m„<^2«2/2) 



n=l 



00 



00 



(All) 



(A12) 



where we adopt changes of variables: y = x — gjuJl and = g„ — c„0(x)/m„w^. Equations 
(lATl) . (!A9l) and (!AT2!1 lead to 

" 27r 



which is independent of a, 6 and r. From Eqs. ( 165|) and ( lAlSp . we finally obtain 



AF = -4ln^^^(^) 



/3 V ^5(0)7 2a;2 



(A13) 



(A14) 
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which is given by Eq. ( 167|) . Equation f lA14p is the same as that for a hnear couphng of 
(f){x) = X in Ref. |47|. 
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FIG. 1: (Color online) Marginal PDFs of (a) (b) P{p), (c) P{u) and (d) P(C) for 

{a,b,T) = (1,0,10.0) (solid curves), (1,0,0.1) (dashed curves), (0,1,10.0) (chain curves), and 
(0, 1,0.1) (dotted curves) calculated by simulations with ajg = 1.0, 70 = 1.0 and ksT = 1.0. P{x) 
and P{p) in (a) and (b) are indistinguishable among the four sets of {a,b,T). P{u) and P{C) for 
{a,b,T) = (1,0,0.1) and (0,1,0.1) are multiplied by a factor of five in (c) and (d) where PDFs are 
arbitrarily shifted for a clarity of figures. 

FIG. 2: (Color online) Time courses of fj,x{t) for (a) r = 0.1 and (b) r = 10.0 with (a, 6) = (1,0) 
and (0, 1) when sinusoidal forces f{t) with g = 1.0 and Tq = 10.0 shown by bottom curves are 
applied. Results for (a, 6) = (0,1) and (1,0) are shifted by five and ten, respectively, for a clarity 
of figures. 

FIG. 4: Lissajous plots of Hx{t) versus f{t) for (a) Tq = 10.0, (b) Tq = 5.0, (c) Tq = 4.0 and (d) 
To = 3.0 with (a, 6, r) = (1, 0, 10.0). 



FIG. 5: (Color online) Time courses of fluctuation of Pxit) for sinusoidal forces with Tq = 6.0, 
10.0 and 100.0 for (a) r = 0.1 and (b) r = 10.0 with multiplicative noise of (a, 6) = (1,0), results 
being successively shifted by two for a clarity of figures. 



FIG. 3: (Color online) The ujq dependence of the output magnitude of Ixi^^o) for sinusoidal forces 
with (a) r = 10 and (b) r = 0.1 obtained by simulations (Sim.) for additive noise of {a,b) = (0, 1) 
(chain curves) and multiplicative noise of {a,b) = (1,0) (solid curves), dashed curves denoting 
theoretical (Th.) results for additive noise [Eq. ([55]) ]. Note that ordinates of (a) and (b) are in the 
logarithmic and normal scales, respectively. 
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FIG. 6: (Color online) The ojq dependence of the stationary fluctuation of pxs for sinusoidal 
forces with (a) r = 10.0 and (b) r = 0.1 for multiplicative noise of (a, 6) = (1,0) (solid curves) 
and additive noise of (a, 6) = (0,1) (dashed curves), error bars expressing variations in p^s- For 
additive noise, pxs = 1-0 independently of loq and r. 

FIG. 7: (Color online) The r dependence of the stationary fluctuation of pxs for sinusoidal forces 
with To = 5 (solid curve) and Tq = 10.0 (chain curve) for multiplicative noise of (a, b) = (1, 0) and 
additive noise of (a, 6) = (0, 1) (dashed curve), error bars expressing variations in p^s- 

FIG. 8: Time courses of Px{t) for (a) additive noise of (a, b) = (0, 1) and (b) multiplicative noise 
of (a, 6) = (1,0) with r = 0.1, 1.0 and 10.0, an applied step force being plotted at bottoms. 

FIG. 9: Time courses of Hx{t) for four sets of {a,b) = (0,1), (0.2,0.8), (0.5,0.5) and (1,0) with 
T = 10.0, an applied step force being plotted at the bottom. 

FIG. 10: (Color online) WDFs of P{W) for ramp forces of (a) r/ = 100.0, (b) r/ = 10.0 and 
(c) Tf = 1.0 with (a,6,r) = (1,0,10.0) (solid curves), (1,0,0.1) (dashed curves), (0,1,10.0) (chain 

curves), and (0, 1,0.1) (dotted curves). Four WDFs in (c) are indistinguishable. 

FIG. 11: (Color online) The r/ dependence of (a) pw, (b) crw, (c) rjw and (d) R for {a,b,T) = 
(1, 0, 10.0) (solid curves), (1, 0, 0.1) (dashed curves), (0, 1, 10.0) (chain curves) and (0, 1, 0.1) (dotted 
curves). Marks denote simulation results and curves are plotted only for a guide of the eye. The 
arrow along the right-hand side ordinate in (d) expresses AF (= —2.0): R = AF when the JE 
holds. 
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FIG. 12: (Color online) WDFs of P{W) with ksT = 1.0 (solid curve), 5.0 (dashed curve) and 
10.0 (chain curve) for a ramp force with r/ = 10.0 for multiplicative noise of (a, b, r) = (1, 0, 10.0). 
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